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ABSTRACT 

PSR J102347.6+003841 is a radio pulsar system with a spin period of 1.69 ms and an 
orbital period of 4.75 hours. Uniquely, it undergoes periods of transient accretion from 
its companion star: it occupies an important position in the evolutionary track from 
X-ray binary to isolated millisecond radio pulsar. Here we present a spectroscopic 
study of this system showing late-type absorption features which match those of a 
G2V star. We find a semiamplitude of 286 ± 3 kms“^ and a best fit orbital period of 
0.1980966(1) days. We combine these measurements with optical photometry which 
suggests the secondary star may be underfilling its Roche lobe by between 15% and 
20%. We weakly constrain the mass of the neutron star to be ^ 2.2 Mq at the 2a level. 
We also discuss the possible origins of the Ha emission line in our template subtracted, 
averaged spectrum. Finally we present and discuss new optical photometry of J1023 
taken during the recent outburst of the system. 

Key words: stars: binaries (including multiple): close - stars: pulsars: individual: 
PSR J102347.6+003841 


1 INTRODUCTION 

Low mass X-ray binaries (LMXBs) are a class of inter¬ 
acting binary star in which material accretes onto a neu¬ 
tron star or black hole from a low mass companion star. 
It is believed millisecond pulsars (MSP) are formed as a 
result of this accretion process spinning up the neutron 
star. The discovery o f accreting millisecond X-ray pulsars 
dWiinands et al. 1998tl confirmed this, however a key ques¬ 
tion is how do these systems evolve into isolated radio mil¬ 
lisecond pulsars? 

PSR J102347.6+003841 (hereafter J1023) is a radio pul¬ 
sar with a spin period of 1.69 ms and an or b ital p eriod of 4.75 
hours. It was discovered by I Bond et al.l il2002ll as part of 
the “Faint Images of the Radio Sky at Twenty Centimeters" 
(FIRST) survey. The spectrum was dominated by a blue 
continuum, much like that of a Cataclysmic Variable (CV) 
star in quiescence, but with stronger than normal Hel and 
Hell lines, suggesting J1023 may be a magnetic CV. They 
also realised this system lay in the region of sky covered by 
the Sloan Digital Sky Survey (SDSS). They analysed the 
archival SDSS data and found the emission lines were dou- 
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ble peaked, indicat ive of the presence of an accretion disc. 
IWoudt et al.l ll2004l l found a repetitive modulation of V ~ 
0.45 mag about a mean of V ~ 17.5 and an orbital period 
of 4.75 hours. They also detected low amplitude flickering, 

suggesting on-going mass transfer, _ 

_ J1023 was observed again bv lThorstensen fc Armstrong! 

(l2005ll (hereafter TA05) where they found only late-type ab¬ 
sorption features, suggesting a dramatic change in the sys¬ 
tem. They combined radial velocity estimates and optical 
photometry which lead them to conclude it was likely that 
the primary star had a mass greater than the Chandrasekhar 
limit and so it was actually a neutron star or black hole, 
rather than a white dwarf. This hypot hesis was supported 
by the optical and X-ray observations of lHomer et al.l ll2006l l 
of J1023, also taken while in quiescence, as they were unable 
to explain the observed X-ray emission in any other way. 

In 2009 Archibald et al. detected radio pulsations from 
the system which confirmed the primary was indeed a neu¬ 
tron star and as such J1023 was confirmed as a LMXB. They 
combined their binary parameters with the radial velocity 
measurements of TA05 to estimate the mass ratio of the sys¬ 
tem, q, to be 7.1 + 0.1. Radio eclipses were also detected, and 
as the line of sight between the pulsar and Earth does not 
intersect the Roche lobe of the companion, they concluded 
these must be caused by material being forced off the com- 
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Figure 1. The average of all 26 spectra is shown below and the 
G2V template used for cross-correlation is shown on top. The 
two most prominent absorption features, the Ca+Fe blend and 
H-alpha lines, are labeled. The data have been shifted to the rest 
frame of the secondary. 

panion star by a pulsar wind, as seen in other similar systems 
llD’Amlco et al.l il200ll l. [Hessels et al.l (l201l[l l. 

In June 201 3 J1023 became radio quiet 

llStappers et al. 2013f l. suggesting, once again, that the 
system had undergone a transition in state. They saw a 
dramatic increase in the 7 -ray flux which along with the 
softer X-ray spectrum of IPatruno et al.l ll2014ll and the 
reemergence of dou ble peaked emission lines as seen by 
iHalpern et al.l (l2013l l confirms accretion, specifically via a 
disc, has begun again in J1023. 

The discovery of J1023 as a millisecond pulsar when 
in quiescence makes the system extremely interesting, es¬ 
pecially given the double peaked emission lines seen in 2002 
and from last year. The history of the system suggests J1023 
undergoes periods of transient accretion from its companion 
star and as such it occupies an important position on the 
track from LMXB to isolated MSP. Here we present analysis 
of archival optical data, combined with new optical photom¬ 
etry, and attempt to more tightly constrain the dynamical 
properties in this important system. 


2 OBSERVATIONS AND DATA REDUCTION 

The Isaac Newton Group (ING) archival data of Rodriguez- 
Gil and Santander presented here were taken using the Inter¬ 
mediate dispersion Spectrograph and Imaging System (ISIS) 
optical spectrometer on the 4.2m William Herschel Tele¬ 
scope at La Palma, Ganary Islands. The red arm of ISIS 
is equipped with a 4096 x 2048 pixel array with an image 
scale of 0.22 arcsec pixeP^ and the high resolution R1200R 
grating with a dispersion of 0.26 A/pixel was used. 

The observations consisted of 26 600s individual expo¬ 
sures obtained during the night of April 12 2009 and cover 
the whole orbital period. We combine this with the pho¬ 
tometry of TAOS to better constrain the binary parameters. 
Data were reduced using the standard IRAljj routines for 

^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 


Name 

Spectral Type 


HR6538 

GIV 

34.023805 

BD+08 2015 

G2V 

21.4970827 

HR1262 

G5V 

37.3525085 

HR159 

G8V 

38.4522552 


Table 1. The range of G-type stars used in the optimal subtrac¬ 
tion process described in the main text. We note the minimum 
occurs for a G2V star. 

bias and sky subtraction, flat fielding and wavelength cali¬ 
bration. The data were then exported to the MOLL10 package 
where they were rebinned onto a common velocity scale for 
analysis. 

We also obtained photometry of J1023 in the V-band 
using the 1.8m Vatican Advanced Technology Telescope 
(VATT) located at Mount Graham International Observa¬ 
tory (MGIO) on April 29 2014. The images were taken using 
the VATT4K GCD with a typical exposure time of 15s. A 
total of 329 images were taken. The images were reduced 
using the standard IRAF routines for flat fielding and bias 
subtraction, and the relatively uncrowded field was analysed 
with the PHOT command in the DIGIPHOT package. We also 
obtained three nights of V-band photometry in June 2014 
using the STE3 GCD on the 1.0m telescope at the South 
African Astronomical Observatory (SAAO). A total of 63 
300s exposures were taken and reduced and analysed in the 
same way as the VATT data. In both cases differential mag¬ 
nitudes were found using two stars, from TAOS, in the same 
field of view as the pulsar: a V = 14.86 magnitude compar¬ 
ison star and a V = 17.40 magnitude star as a check. 


3 RESULTS 

Figure [T] shows the average pulsar spectrum along with our 
choice of template. The spectral type of the secondary star 
was estimated by performing a fit to our selection of G- 
type stars. The results are displayed in Table [T] 

The radial velocities of the individual spectra were 
found through cross-correlation with a high signal-to-noise 
G2V template spectrum observed on the same night as the 
target spectra. The wavelength range used was 6200A - 
68OOA (excluding the Ha region) and the cross-correlation 
was carried out using the xcor routine in molly and fxcor 
in IRAF. Figure shows the radial velocity curve obtained 
through the cross-correlation technique. A sine wave fit 
yields a systemic velocity of 7 = 0.45 ± 1.6 kms“^, in good 
agreement with TA05, and a semi-amplitude of 286 ± 3 
kms”’^ for K 2 compared to 268 ± 4 kms“^ found by TA05. 
However our value of K 2 supersedes theirs as our spectra are 
of a higher resolution (lA compared to their 3.5A) and due 
to the better sampling of our radial velocity curve. 

In order to refine the or bital period of TA05 we used th e 
O-C method described by iBevington fc RobinsonI ll2003ll . 
This resulted in a best fit orbital period P = 0.1980966(1) 

Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation 

^ http://www2.Warwick.ac.uk/fac/sci/physics/research/astro/ 
people/marsh/software 
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Figure 2. Radial velocity curve folded over the best fit period 
of 0.1980966(1) days. Error bars have been scaled to achieve a fit 
with a reduced = 1- 


Free Parameters 

Results 

/ 

U.OO_0.02 

i 

54°+i: 

Te// 

5655 ± 85 K 

La; 

1033 erg s-i 


Table 2. Results of the ELC modelling for the secondary star 
underfilling its Roche Lobe. 


days at the Icr confidence limit. This improves on the 
value found by TAPS, however not on the value found by 
[Archibald et al.l ll2009ll of 0.1980962019(6) days which was 
obtained through timing measurements of radio pulsations 
of the pulsar. 

Using our K 2 and orbital period we estimate the mass 
function, f(m 2 ), to be 0.48 ±0.02 Mq. 



Figure 3. Results of our Roche lobe underfilling modelling of 
J1023 using ELC and the V-band photometry of TA05. Data are 
repeated for clarity. 


4 AN ECLIPSING LIGHT CURVE MODEL 

We next attempted to constrain the system parameters 
by utilising the E clipsing Light Curve Code, (ELC), of 
lOrosz et al. (l2000l ). to simultaneously fit the V-band light 
curve of TA05 and o ur new radial velocity curve. A Monte 
Carlo Markov Chain llTeemark et al. 20041 ) was used to de¬ 
termine errors. The bounded free parameters were the Roche 
lobe filling factor, f, the effective temperature of the sec¬ 
ondary, 5570K < Te// < 5750K (in line with a G2V spec¬ 
tral type), the irradiating luminosity of the pnlsar, 10 ^^ - 
4 X10®^erg s“^ (where the upper limit is determined by 
the spin down luminosity of the pulsar), the mass ratio of 
the system, 7.0 < q < 7.2 (we adopted the mass ratio of 
[Archibald et al.l ([2009[ ). 7.1 ± 0.1), and the inclination of 
the system, 36° < i < 55°. The upper b ound on the in¬ 
clinat ion comes from previous estimates of [Archibald et ^ 
|[2009l) while the lower bound comes from noting that, for 
a fixed K 2 , period and mass ratio, the calculated mass of 
the primary is entirely dependent on the inclination of the 
system, and an inclination of less than 36° gives a primary 
mass greater th an 3 M( 7 ), the maxim um theoretical mass of 
a neutron star |[Chamel et al. 2013[) . The albedo of the sec¬ 
ondary star was set to 0.5, in keeping with the star being a 
low mass (e.g. --^0.2 Mq for a 1.4 Mq pr i mary ), convection 
dominated star as noted bv [Deller et al.[ ([2012[ ). 

The results can be seen in Table with the fit to the 
V-band shown in Figure[31 The left panel in Figure [T] shows 
the correlation between the inclination of the system and 
the filling factor of the companion. Th e inclination here i s 
higher than the 44° ± 2° inferred from [Peller et al.l (|2012[ ) 
however they have a higher mass for the neutron star which 
assumed a Roche lobe filling secondary. As we have allowed 
the Roche lobe filling factor of the secondary star to vary, 
we are only able to provide an upper limit to the mass of 
the neutron star. The limits are 1.4 Mq at the Icr level and 
2.2 Mq at the 2cr level. This can be seen in the right panel 
in Figure [4] which shows the correlation between the mass 
of the neutron star and the filling factor of the companion. 

It is worth noting that our modelling assumes the sec¬ 
ondary star has a convective envelope which, given the evo¬ 
lutionary history of the system, may not necessarily be the 
case. We note that a radiative envelope changes the optical 
modulation by approximately 0.1 magnitudes which would 
have a significant effect on the parameters estimated here. 

We also investigated the effects of veiling on the light 
curve (for example, due to a constant contribution from an 
accretion disk). For a Roche lobe filling factor of one we 
require a relatively high veiling (of order 100 %, or twice the 
observed minimum flux) to best match the observed data. 
However, this is not consistent with veiling constraints we 
obtain from spectroscopic measurements: these provide an 
upper limit of 30% for a G2V template. 

By using a temperature of 5655K for a G2V star and 
radius of the companion star as computed by ELC of 0.39 ± 
0.07 Rq, we found the absolute magnitude of the companion 
should be 6.88 ± 0.07, which is in good agreement with the 
absolute magnitude of 6.84T0.07, c alculated using a dis tance 
to the system of 13681(39 Pi- from [Peller et al.l l[ 2012 l ) and 
using an extinction of 0.14 in the V band. 
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Figure 4. The left panel shows the inclination of the system against the Roche lobe filling factor of the secondary star and the 
right panel shows the mass of the primary against the filling factor. The contours here represent constant values of reduced 
and the dashed lines are the Icr confidence limits. The minimum occurs at 135.88 which corresponds to a reduced x^=l- 
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Figure 5. The residual Ha emission seen as a result of the opti¬ 
mal subtraction. 


4.1 Residual Ha Emission 

The cross correlation, averaging of the J1023 spectra and 
optimal subtraction of the template spectrum resulted in a 
residual Ho emission line with equivalent width (—1.3±0.2) 
A(where the error was found by varying the continuum level 
when fitting the emission line) and can be seen in Figure 
[5l To estimate the K 2 of the Ha emission line we cross 
correlated and optimally subtracted the template spectrum 
from each of our 26 individual exposures to produce a ra¬ 
dial velocity curve of the emission line. We found a K 2 of 
290 ± 4 kms“^ which is consistent with the absorption line 
value (286 ± 3 kms“^), which suggests that the source of 
the residual emission Is from the secondary Itself, possibly 
an irradiated inner hemisphe re. Alternatively as discussed in 
iHernandez fc Casared (l2010l l. and the references therein, the 
measured equivalent width is consistent with Ha emission 
due to a chromospherically active secondary star. 


5 A CHANGE IN STATE 

As stated in Section [U J1023 reverted back to its accre¬ 
tion phase in June 2013. We obtained optical photometry 
of the system over three nights using the Im telescope at 
the South Africa Astronomical Observatory (SAAO). This 
was combined with a single night’s observation at the Vati¬ 
can Advanced Technology Telescope (VATT) at the Mount 
Graham International Observatory (MGIO) in Arizona. To¬ 
gether these observations cover three quarters of the or- 
bital period and we combine this with the V-band data of 
iHalpern et al.l (l2013l l which gives us complete phase cover¬ 
age (Figure [ 6 l) 

The data indicate that the system has brightened by 
over a magnitude compared to the average V-band magni¬ 
tude of ~17.5 when in quiescence. 

The data suggest two types of behaviour in this state: 
one appears as a somewhat smooth, albeit asymmetric light 
curve, with minimum at phase zero, whereas during the 
epoch of the VATT observations, a second type of more 
rapid, flaring variability is observed. 


5.1 Comparison to Similar Systems 


5.1.1 Long Term Variability 

XSS J12270-4859 (hereafter J12270) was discovered as part 
of the Rossi X-ray Timing Explorer Project (RXTE) slew 
survey |Sazona^_et_^l_2013 ) and was initially classified as a 
GV bv iMasetti et al.l 1 200a) due to the presence of optical 
emission lines. Furthe r multi-wavelengt h obse rvations cast 
doubt on this, indeed Ide Martino et aO ll2010ll noticed the 
position coincide d with a Fe r mi 7 -r ay source (IFGL J1227.9- 
4852) which lead iHill et al.l ll201lll . amongst others, to sug¬ 
gest t hat J12270 may ha rbour an active radio MSP. In their 
paper iBassa et al. liiH) describe changes in the system as 
it transitions from an accreting LMXB phase to that of a ro¬ 
tation powered MSP. They note that the optical light curve 
decreases by between 1.5 and 2 magnitudes while there is a 

















































Roche Lobe Underfilling in PSR J102347.6+003841? 5 



Figure 6. The three sets of V-band photometry taken while the 
system was in its recent outburst phase. Figure I6al shows the 
SAAO data (black open squares) taken using the Im telescope on 
the nights of 2014/06/06, 2014^6/12 and 2014 /06/13 and Atel 
data (red filled circles) of iHalpern et al] ll2013tl taken using the 
1.3m McGraw-Hill telescope at the MDM Observatory on Kitt 
Peak. Figure [60 shows the VATT data taken on the night of 
2014/04/29. Data are repeated for clarity. 

corresponding decrease, by a factor of 10, in the X-ray count 
rate and a decrease in the 7 -ray brightness, again by a factor 
of between 1.5 and 2. This matches J1023 very well as in its 
change of state from MSP to LMXB the optical luminosity 
increases by between 1 and 2 magnitudes from quiescence, 
the X-ray luminosity increases from ~10^^ erg s“^ to ~10®^ 
erg s“^ while the 7 -ray flux increased by a factor of five. 

5.1.2 Short Term Variability 

Figure I 6 bl shows that some of the flaring events occur on 
timescales shorter than the sampling time of 15s, while 
other, higher amplitude, flares last 100s. While in similar 
systems both the short and l ong timescale variat ions tend 
to last longer than in J1023, iHvnes et al.l ll2002h . in their 
study of V404 Cyg, suggest that shorter timescale flickering 
may be a common feature of LMXB systems while the larger 
amplitude events are more difficult to explain, but may be 
due to instabilities in the accretion disc. Si milar short and 
long t imescale events are seen in Cen X-4 and lShahbaz et aP 
1I2OIOII compared the observed colours to those predicted by 


three emission mechanisms: blackbody, synchrotron radia¬ 
tion and an optically thin layer of Hydrogen. They favour 
this third mechanism and attribute the source of the flares 
to the accretion disc. However the amplitude of the flaring 
we observe in J1023 appears unusual in comparison to other 
blac k hole or neutron star related systems. 

[Archibald et al.l ([20iri l found evidence of X-ray emis- 
sion modulated at the rotational period of the pulsar. 
[Russell et al.[ ([2007['l analysed optical/infrared (OIR) and X- 
ray data from 19 Neutron Star X-ray Binaries (NSXB) to 
compare the observed OIR and X-ray fluxes to their ex¬ 
pected values if the OIR emission is dominated by thermal 
emission from a X-ray heated disc or synchrotron emission 
from the jet. They found that thermal emission due to X- 
ray reprocessing could explain all their data except at high 
luminosities. Synchrotron emission from the jet was found 
to dominate the NIR light curve above La, ~ 10^® erg s“^ 
and the optical above Lj, ~ 10^^ erg s“^ . This suggests that 
the flaring seen in our VATT data is not due to synchrotron 
emission, given the spin down luminosity of the pulsar, Lj, 
~ erg s"^ 


6 CONCLUSIONS 

The results of our modelling suggest that the secondary star 
in J1023 is underfilling its Roche lobe between 15% and 
20%, albeit at the Icr level. If confirmed, this suggests that 
the mode of accretion in this system is more likely to be 
via a wind from the secondary star, rather than Roche lobe 
overflow which is more commonly the case. This may in 
turn be related to the accreting/non-accreting transitions 
that have been observed from this system. 

We also observe evidence for Ho in emission in the qui¬ 
escent state, which we associate with either the irradiation 
of the secondary star by the spin down luminosity of the 
pulsar, or from an active chromosphere on the secondary. 

As our photometric observations of J1023 in this out¬ 
burst phase have been taken over several epochs, in which 
only part of the orbital phase has been sampled, more sys¬ 
temic, multi-wavelength, observations are required. Further 
multi wavelength observations once it returns to quiescence 
are encouraged, especially phase-resolved spectroscopy and 
more accurate photometry to confirm the underfilling of the 
Roche Lobe of the secondary star, the source of the Ha 
emission, and more tightly constrain the mass of the neu¬ 
tron star. 
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